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The Hawking evaporation of small black holes formed by the collapse of dark matter at the center 
of neutron stars plays a key role in loosing the constraint on the mass of asymmetric bosonic non- 
interacting dark matter particles. Different from previous works we show that such a kind of dark 
matter is viable in the mass range from 3.3 GeV to ~ 10 TeV, which covers the most attractive 
regions, including the preferred asymmetric dark matter mass ~ 5.7 GeV as well as the 5—15 GeV 
range favored by DAM A and CoGeNT. 



PACS numbers: 95.35.+d, 95.30Cq, 97.60.Jd 



I. INTRODUCTION 



There is compelling evidence from astrophysical and 
cosmological data that the dominant component of the 
matter in the universe is in the form of "dark matter" 
that only interacts very weakly with ordinary matter [l[ . 
One of the remarkable features of the standard cosmologi- 
cal model is that the dark matter density today is compa- 
rable with that of the baryons (i.e., pdm ~ 5.7pbaryon) H- 
Such a fact could be naturally explained by an asymme- 
try in dark matter similar to that in baryons Q . Though 
higher or lower masses are possible [3, Q , the asymmet- 
ric dark matter (ADM) models prefer to have DM mass 
~ /ODMWn/yObaryon ~ 5.7 GeV [6], which is within the 
inferred mass range (5 — 15 GeV) of the dark matter can- 
didates found by DAMA Q and CoGeNT H, where m„ 
is the mass of neutrons. In the ADM scenarios, despite 
the DM particles can in principle have weak-scale inter- 
actions and therefore sizable scattering cross sections off 
baryons, annihilations do not occur because of the pres- 
ence of an asymmetry in the DM sector between particles 
and anti-particles. Therefore the ADM particles can be 
captured by astrophysical bodies such as the sun, earth, 
white dwarfs and neutron stars 0-[H| ■ The captured par- 
ticles will eventually be thermalized and centered in the 
core of the star, namely, a dense dark matter core forms 
[9l-ITTj . Due to the non-annihilation of ADM particles, 
the build up of a concentration at the center of neutron 
stars may lead up to a collapse of these particles into a 
small black hole and the neutron star may be destroyed. 
If the central temperature of very old isolated neutron 
stars falls below the critical temperature to form a Bose- 
Einstein condensate (BEC), the bosonic particles in the 
ground state condense and no longer follow the thermal 
distribution. Consequently the process of collapse will 
be significantly enhanced. Therefore the existence of very 
old neutron stars in the Galaxy imposes a valid constraint 
on non-annihilating bosonic weakly interacting massive 
particles (WIMPs). In the initial work by Goldman and 
Nussinov [lfj a mass larger than ~ 1 GeV has been ex- 



cluded. In other words, the preferred ADM models with a 
mass ~ 5.7 GeV can not be non-annihilating and bosonic. 
Such a result is a piece of evidence that the compact ob- 
jects can be valuable probe of the physical properties of 
dark matter particles. Recently, after taking into account 
the Hawking evaporation of the newly-formed small black 
hole [H, UH , Kouvaris and Tinyakov [3] have shown that 
the non-interacting ADM particles can be bosonic as long 
as they have a mass m x > 16 GeV, i.e., the allowed mass 
range of non-interacting bosonic ADM particles has been 
considerably extended. Such a progress is remarkable but 
the 5 - 15 GeV range favored by DAMA and CoGeNT has 
still been excluded. In this work we show that with the 
advanced estimate of Hawking evaporation rate [l3j as 
well as the rate of accretion onto the black holes formed 
at the center of neutron stars, asymmetric bosonic non- 
interacting dark matter with a mass > 3.3 GeV is viable, 
which covers the most interesting regions, including the 
preferred ADM mass ~ 5.7 GeV as well as the 5—15 GeV 
range favored by DAMA and CoGeNT. We also find that 
the asymmetric bosonic non-interacting dark matter par- 
ticles with a mass above ~ 10 TeV has been excluded by 
current neutron star observations. 



II. CONSTRAINT ON THE MASS OF 
ASYMMETRIC NON-INTERACTING BOSONIC 
DARK MATTER 

WIMPs scatter with nuclei and then lose their kinetic 
energy when pass through a star. If the kinetic energy is 
smaller than the stellar gravitational potential the dark 
matter particles will get captured [9( and concentrate in 
the center within the radius [151 



rth 



74 cmTjf 



lOOGeV 



U/2„-l/2 
I Pc.15 i 



(1) 



where T c (p c ) is the temperature (density) of the 
star core. Through this work we denote p c ,i5 = 
p c /10 15 g cm" 3 and T c , 6 = T c /10 6 K. 
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The mass of the ordinary matter contained in such a 
volume is 



and 



then 



nwf « 10 4a GeV T c 2 6 ( 
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The mass of the captured dark matter can be estimated 
as [HI 

m dm « 1.3 x 10 39 GeV /( /*" )(^_), ( 3 ) 

0.3 GeV cm J 10° yr 

where pdm is the energy density of dark matter surround- 
ing the neutron star, / = min{l, o x -^.-j\ and o~ x is the 
spin-independent cross section of the dark matter scat- 
tering with the neutron/proton. 



A. Excluding the light asymmetric bosonic dark 
matter 

For bosonic WIMPs, a Bose-Einstein condensate 
(BEC) may be formed. Since such a state is more 
compact, the self-gravitation regime starts much earlier. 
The particle density required to form BEC is tibec — 
1.5 x 10 33 cm' 3 ( 10 q^ cV ) 3/2 ( ) 3/2 and the total mass 
of the needed particles is 



2.5 x 10 41 GeV ( 



100 GeV 



-3/2 



(4) 



The size of the condensed state is determined by the 
radius of the wave function of the WIMP ground state in 
the gravitational potential of the star, i.e., 

r cond «1.6xlO- 5 cm(H^X) 1 /2 ; (5 ) 



which is much smaller than rth- The self- gravitation 
regime reaches for a mass of the condensed matter 



m-BEC,sclf 



: 8 x 10 24 GeV (m x /100 GeV)- 3/2 p c ,i5- 



For non-interacting bosons only the uncertainty prin- 
ciple opposes the collapse. As long as the mass of self- 
gravitating bosons reaches flij 
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they will collapse into a black hole, where m p i = 1.22 x 
10 19 GeV is the Planck mass. For m x > 1 keV that is of 
our interest, m c 3> ^BEC.seif • Hence it is m c rather than 
TOBEC.seif that provides a key constraint on the mass of 
bosonic non-interacting ADM. 

For a black hole to form, the condition either 
m c < m B EC < "idm or m se if < m dm should be 
satisfied. The request tobec < Wdm suggests that 
the neutron star should have an age iNS > 2 x 

10 8 ( I5 ^ 5V )I?, 6 p-; 5 /2 /- 1 ( o.3 g^-^ )- 1 - While 
m c < mBEC requires 



m x > m\o 
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10 5 yrs T% 2 p:% A f- X 3G PJ? cm _» )-\ Current observa- 
tions of nearby highly magnetized neutron stars suggest 
that these objects usually cool down to a surface tem- 
perature T s ~ 10 6 K (correspondingly T c ~ 10 8 K if the 
relation T c w 10 8 K [(T s /10 6 ) 4 /(5/2 x 10 14 cm s" 2 )] 455 
[l8T | still holds, where g is the surface gravity of 
the neutron star) in a timescale t ~ 10 6 yrs [l9l |. 
If such objects can be identified in regions with 
Pdm > 120/- 1 GeV cm" 3 , one has m low = 0.2 MeV. 

For ADM particles lighter than ~ 0.2 GeV, it 
is rather hard to have m se if < mdm- For ex- 
ample even with very extreme parameters pdm ~ 
3 TeV cm" 3 , t - 10 10 yrs and T c - 10 4 K, we have 
mdm ~ 10 47 GeV, which is still smaller than m se if ~ 
10 48 GeV (T c /10 4 K) 3 / 2 (0.1 GeV/m x ) 3 / 2 . Therefore cur- 
rent neutron star observations are not against the keV to 
sub-GeV non-interacting bosonic ADM and the future 
observations of highly magnetized neutron star in dense 
dark matter region may exclude the mass range from 
sub-MeV to sub-GeV. 

After its formation, whether the black hole can sur- 
vive or not depends on the competition between Hawking 
evaporation and simultaneous accretion of the surround- 
ing material. Let's discuss these two processes in some 
detail. 

The Hawking radiation mimics thermal emission 
from a blackbody with a finite size and a temperature of 
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where Mbh is the mass of the black hole. The mass-loss 
rate of an evaporating black hole can be expressed as 



dMBH = -5.34 x 10 25 g s- 1 F(Af BH )( MBH ' 1 



dt 
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(9) 



where F(Mbh) accounts for degrees of freedom of each 
species of radiated particles and can be estimated with 
eq.(9) of [llj]. When the temperature of the black hole ex- 
ceeds the quark-gluon deconfinement temperature (~ 200 
MeV), the quark, gluon and lighter particles will be emit- 
ted. Hence the mass loss rate can be considerably larger 
than the estimate dM B }i/dt = -& 4 /153607rG 2 M| H 
adopted in both [TJ[ and [lj|, where G is Newton's con- 
stant, c is the speed of light and h is the reduced Planck's 
constant. 

Bondi accretion onto black hole has been widely 
discussed since 1952 [17[. For current purpose we need 
to take the relativistic approach. Following Appendix 
G of ,2(|, the rate of accretion of rest mass onto the 
black hole is Airm n nur 2 — M = constant (indepen- 
dent of the radius r to the black hole), where n is the 
baryon number density and u is the inward radial ve- 
locity component. To calculate an explicit value for M, 
we need to adopt an equation of state and take it to be 
P = Kn T , where P is the pressure. In K and T 
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are assumed to be constant. Currently the density is so 
high that the equation of state likely changes when n 
reaches a critical value. For n < n t — 1.15 baryon fm -3 
(pt w 1.9 X 10 15 g cm -3 ), motivated by the latest as- 
trophysical measurement on the equation of state of 
the ultra-dense material in neutron stars (2l| we take 
P k, 2 MeV fm~ 3 (n/0.2 baryon fm~ 3 ) 3 . At n = n t , 
P w nm n c 2 /3 and quark matter likely forms (the tem- 
perature is denoted as T t ), for which T w 4/3 and the 
sound speed is a s ss c/V3. The energy density (baryon 
number density) is proportional to T 4 (T 3 ). Adopting 
eq.(G.17) of [2Q|, at the critical radius r s rj 3GM Bii /c 2 , 
we have u s = c/y/E, and u 2 /a 2 — 1/2. Hence the rela- 
tivists Bernoulli equation (i.e., eq.(G.22) of [2(|) gives 



Ts /T t = (a 2 Ju 2 a f/ 2 = y/2, n s = (a 2 Ju 2 s ) 3 / 2 n t 
and then 

M = Airm n n s u s r 2 « 24v / 37r / 9 t G 2 M| H /c 3 . 
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Our accretion rate is lower than the preliminary estimate 
with classical Bondi accretion formula [l4| by a factor of 
- 10. 

The black hole will disappear as long as M + dA ^f H < 0, 
such a request defines a critical mass of the black hole, 
i.e., 



A/ B H,c = 5.2xl0 13 g( 



Pt )-l/4 ^( M BH) ^i/4 



2 x 10 15 g cm 
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The Hawking temperature is Tr ~ 203 MeV, for which 
F(Mbh) « 6. 

The mass of the dark matter particles should satisfy 
(i.e., eq.©) 



i x > m x , c = 9.5 x 10 3 ' GeV / (Mbh,cC ) 



3.3 GeV ( 



Pt ^i/ 4 ^(^bh)^i i 



2 x 10 15 g cm 



G 



ttf) 



otherwise the black hole will be bigger and bigger and 
finally the whole neutron star is destroyed. 



B. Excluding the TeV asymmetric bosonic dark 
matter 



Very old neutron stars will be destroyed too if BEC 
does not form (i.e., Mdm < M BEC ) but the self- 
gravitating dark matter core is massive enough for black 
hole to form (i.e., Md m > M Bel! and Md m > 2m 2 j/ '^m x ) . 
The existence of old neutron stars in the Galaxy thus im- 
poses a tight constraint on the heavy ADM model [lj|. 
In this case the request Md m > 2m 2 >l /irm x usually does 
not provide us new information and we focus on the first 
two requests, i.e., 



> mi = 0.5 GeV / 
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and 

m x > m 2 



840 TeV T c>6 /" 2/3 ( 
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Therefore we need m x > m up 



106 yr ; 

max{rni, m-i\. Due to 
the relatively fixed p c and the reverse sign of the indices 
of mi and ni2 on /, /5dm, t, it is straightforward to show 
that m up reaches the minimum when mi = m 2 , for which 
we have 

Pdm \5/12 



T c , crit = 2.8 x 10 4 K f' 12 { 
t 



( 



10 6 yr 



0.3 GeV cm- 3 ' 
^VYf, (15) 



and then 

m up = 24TeV /"*( 



0.3 GeV cm" 3 



' V 10 6 yr 



r 1/8 
'Pc,15- 

(16) 

Since / > 1, for the extreme parameters pdm ~ 
3 TeV cm -3 and t ~ 10 10 yrs, the above constraint yields 
that m up > 240 GeV. However, such a tight constraint 
is not reachable because the required T c is high up to 
~ 10 8 K, which is impossible for a neutron star as old 
as the Universe. Nevertheless, even in such an extreme 
scenario, the resulting m up ~ 240 GeV is far larger than 
m Xi c £3 3.3 GeV, demonstrating that the window for non- 
interacting bosonic asymmetric dark matter models does 
open. 

For isolated neutron stars, with the relation between 
T c and TL, we have 



T s , crit « 1.1 x 10 4 K /°' 23 ( 



Pdr. 



0.3 GeV cm 



— ) c 

-3 I 



(- 



t 



0.23^0.25„l/4 
Pc,15yi4.3- 



(17) 



10 6 yr 

To form a black hole, m x should satisfy m x > 
TOup^/T^crit) -5 - 46 for T s < T S:Crit and m x > 
TOup^/r^crit) 1 - 82 for T s > r S:Crit . Among the param- 
eters involved in the above equation, T s , pdm and t are 
measurable, and p c and g are relatively well determined. 
So the tightest constraint on m x is achievable. 

The cooling of neutron stars is still to be better un- 
derstood. The fit to the surface temperature of the best 
studied isolated neutron stars suggests that for a mass 
- 1.4M e , one likely has T s ~5x 10 5 K (~ 2.3 x 10 4 K) 
at t - 10 6 yrs (~ 2.5 x 10 7 yrs) [II]. With ea.([T7]l we 
have T s , crit w 2.3 x 10 4 K - T s at a time t ~ 2.5 x 10 7 
yrs. The existence of isolated neutron stars older than 
2.5 x 10 7 yrs thus requires that 



m , 



<m up «10TeV 



^0.3 GeV cm- 3 ; ' (18) 
otherwise such stars will have been destroyed by the black 
hole formed at their center (the formed black hole has a 
mass 3> 5.6 x 10 38 GeV, for which the mass loss due to 
Hawking evaporation is too small to balance the Bondi 
accretion. It is straightforward to show this with the 
formulae presented in Section II A). 
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FIG. 1: The exclusion regions as a function of m x and 
<7 X for isolated neutron stars at local DM density pdm = 
0.3 GeV cm -3 . The left exclusion region is the case of J0108- 
1431 (t N s ~ 1-6 x 10 s yrs and T c ~ 10 6 K [H) while the right 
exclusion region is based on eq. lfTH)) . 



III. CONCLUSION 



Bondi accretion rate has been overestimated (see sec- 
tion II A for the details). After the corrections, the 
asymmetric bosonic non-interacting dark matter with a 
mass > m XjC « 3.3 GeV is found to be viable. We 
also show that for m x > m up ~ 10(pdm/0.3 GeV cm~ 3 ) 
TeV, the black hole can form (without reaching the Bosc- 
Einstein condensation regime) and swallow the neutron 
star. Interestingly, even in the most-extreme/unrealistic 
scenario, we have m up ~ 240 GeV, which is much larger 
than m XjC ~ 3.3 GeV, demonstrating that there is indeed 
a nice window for non-interacting bosonic asymmetric 
dark matter models (see section II B for the details). We 
conclude that the asymmetric non-annihilating bosonic 
dark matter is viable if its mass is in the range from 3.3 
GeV to ~ 10 TeV, which covers the most attractive re- 
gions, including the preferred asymmetric dark matter 
mass ~ 5.7 GeV as well as the 5—15 GeV range fa- 
vored by DAM A and CoGeNT (see FigHJ. Moreover, 
we find out that current neutron star observations are 
not against the keV to sub-GeV non-interacting bosonic 
ADM and the future observations of magnetized neutron 
star in dense dark matter region may exclude the mass 
range from sub-MeV to sub-GeV (see eq.© and the sub- 
sequent two paragraphs for the discussion). 



For the bosonic non-interacting asymmetric dark mat- 
ter particles, the build up of a concentration at the cen- 
ter of neutron stars may lead up to a collapse of these 
particles into a small black hole and the neutron star 
may be destroyed. Therefore the existence of old neu- 
tron stars in the Galaxy imposes a tight constraint on 
such a kind of dark matter, as firstly noticed in 1989 
10]. A robust constraint is however not achievable until 
the Hawking evaporation of the newly-formed small black 
hole and the simultaneous accretion have been properly 
addressed. We find out that in previous works the Hawk- 
ing evaporation rate has been underestimated while the 
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